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Abstract 

! Transverse emission of preequilibrium nucleons, light clusters (complex particles) and charged 
pions from the isotopic ^^^'^^^Sn+^^^'^^^Sn reactions at a beam energy of 400^ MeV, to extract the 
high-density behavior of nuclear symmetry energy, are investigated within an isospin and momen- 
tum dependent transport model. Specifically, the double ratios of neutron/proton, triton/helium-3 

■ and tt^/tt^ in the squeeze-out domain are analyzed systematically, which have the advantage of 
reducing the influence of the Coulomb force and less systematic errors. It is found that the 
transverse momentum distribution of isospin ratios strongly depend on the stiffness of nuclear 

] symmetry energy, which would be a nice observable to extract the high-density symmetry energy. 

; The collision centrality and the mass splitting of neutron and proton in nuclear medium play a 
significant role on the distribution structure of the ratios, but does not change the influence of 
symmetry energy on the spectrum. 
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Heavy-ion collisions with neutron-rich beams at intermediate and relativistic energies are a 

unique tool to extract the information of nuclear equation of state (EoS) of isospin asymmetric 

nuclear matter under extreme conditions, such as high density, high temperature and large isospin 

asymmetry etc. The energy per nucleon in the isospin asymmetric nuclear matter is usually 

expressed as E{p, 6) = E{p, 6 = 0) + EsYmip)^^ + 0(5^) in terms of baryon density p = Pn + Pp, 
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relative neutron excess 6 = {pn — Pp)/{Pn + Pp), energy per nucleon in a symmetric nuclear 
matter E{p,6 = 0) and nuclear symmetry energy -Esym = \s=o- Typically, two different 

forms have been predicted by some microscopical or plienomenological many-body approaches. 
In one, the symmetry energy increases monotonically with the baryon density, and in the other 
the symmetry energy increases initially up to a suprasaturation density and then decreases with 
baryon density. The density dependence of the symmetry energy around and below the normal 
densities has been roughly constrained and investigated from heavy- ion collisions [1] [2] . Moreover, 
the high-density behavior of nuclear symmetry energy predicted by transport models associated 
with the existing experimental data deviates largely because of the inconsistent treatment of the 
mean-field potentials and the unclear in-medium nuclear interaction and properties of resonances 
[21 in El E], and an opposite conclusion was drawn. Furthermore, new experimental data related to 
high-density observables and modifications of transport models including the in-medium effects in 
two-body collisions and in mean-field propagation, are very anticipated. In addition to providing 
reasonable understanding of the reaction dynamics and the structure of neutron-rich nuclides, the 
high-density information of symmetry energy also has an important application in astrophysics 
regarding such subjects as the structure of neutron star, the cooling of protoneutron stars, the 
nucleosynthesis during supernova explosions of massive stars etc |7]. The importance of high- 
density symmetry energy in hadron-quark phase transition was stated and investigated thoroughly 
by Di Toro et al. [8]. 

The preequilibrium neutron and proton transverse emission in isotopic reaction systems at an 
incident energy of 50 MeV/A was measured at the National Superconducting Cyclotron Labora- 
tory (Michigan State University, East Lansing, MI, USA) ^ and analyzed with transport models 
for constraining the symmetry energy at subnuclear densities [lOl [11] , which have been verified 
as a nice probe for extracting the density dependence of the symmetry energy. The double neu- 
tron/proton ratio of free nucleons was investigated in the high-energy range as a probe of symmetry 
energy at high-baryon density using the isospin Boltzmann-Uehling-Uhlenbeck (IBUU04) model, 
which was concluded that the symmetry energy results in a turnover spectrum in comparison to 
the kinetic energy distribution of the lower energy 50 MeV/A [12]. In this work, the transverse 
emission of prompt nucleons and light clusters as well as charged pions produced in the isotopic 
Ii2,i24gj^_l_ii2,i24gj^ reactions, to extract the symmetry energy at supranormal densities, is investi- 
gated within an isospin and momentum dependent transport model (Lanzhou quantum molecular 
dynamics (LQMD)). The influence of symmetry energy and mass splitting of neutron and proton 
on the isospin emission is systematically analyzed. 
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The LQMD model has been updated in order to extract the density dependence of the nuclear 
symmetry energy from heavy-ion collisions. The main inelastic channels of baryon-baryon and 
meson-baryon collisions as well as the decay of resonances for producing pion and strangeness 
particles were implemented in the model [H [T21 [IH [151 [IS]- The momentum dependence of the 
symmetry potential was also included in the model, which results in a splitting of proton and 
neutron effective mass in nuclear medium pTfj. In the LQMD model, the time evolutions of the 
baryons (nucleons and resonances (A(1232), N*(1440), N*(1535))) and mesons in the system 
under the self-consistently generated mean-field are governed by Hamilton's equations of motion, 
which read as 

OH .OH 

Here we only consider the Coulomb interaction for charged hyperons. The Hamiltonian of baryons 
consists of the relativistic energy, the effective interaction potential and the momentum dependent 
interaction. The effective interaction potential is composed of the Coulomb interaction and the 
local interaction 

Uint = Ucoul + Uioc- (2) 

The Coulomb interaction potential is written as 

Ucoui = ^ E ^erf{r,,/VAL) (3) 

where the ej is the charged number including protons and charged resonances. The = |rj — rj| 
is the relative distance of two charged particles. 

The local interaction potential is derived directly from the Skyrme energy-density functional 
and expressed as 

Uioc = J ViMr))dr. (4) 
The local potential energy-density functional reads 



2 po 1 + 7 Po 2po 
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[V{pn-Pp)Y, (5) 



2po 

where the p„, Pp and p = Pn + Pp are the neutron, proton and total densities, respectively, and the 
6 = {pn — Pp)/{Pn + Pp) is the isospin asymmetry. The coefficients gr, Qsur, o'sur are related to the 
Skyrme forces and the parameter Sly6 is taken in the calculation |14j. The bulk parameters a, /3 
and 7 are readjusted after inclusion the momentum term in order to reproduce the compression 



modulus (here, K=230 MeV) and the binding energy of isospin symmetric nuclear matter at 



saturation density, which are taken as -296.6 MeV, 197 MeV and 1.143, respectively. The E[ 



loc 
sym 



is the local part of the symmetry energy, which can be adjusted to mimic predictions of the 
symmetry energy calculated by microscopical or phenomenological many-body theories and has 
two-type forms as follows: 



and 



KUp) = lc,Up/Por% (6) 



KmiP) = asym{p/Po) + bsym{p/Pof- (7) 



The parameters Cgym, dsym and hsym ai'c taken as the values of 52.5 MeV, 43 McV, -16.75 McV 
and 23.52 MeV, 32.41 MeV, -20.65 MeV corresponding to the mass splittings of m* > m* and 
m* < m*, respectively. The values of 7^=0.5, 1., 2. correspond to the soft, hnear and hard 
symmetry energy, respectively, and the Eq. (7) gives a supersoft symmetry energy, which cover 
the largely uncertain of nuclear symmetry energy, particularly at the supra-saturation densities. 
All cases cross at saturation density with the value of 31.5 McV. 

A Skyrme-type momentum-dependent potential is used in the model as follows: 

Umom^ Y~ Yl ^Cr,T'5r,Ti5T',Tj j J J dpdp' dr 

x/,(r, p, 0[ln(e(p - p')' + 1)] Vi(r, p', t). (8) 



Here C^.r — C^omi^ + x), C'r.r' — C^omi^ ~ {'^ '^') ^^d the isospin symbols t(t') represent 
proton or neutron. The sign of x determines different mass splitting of proton and neutron in 
nuclear medium, e.g. positive signs corresponding to the case of m* < m*. The parameters Cmom 
and e was determined by fitting the real part of optical potential as a function of incident energy 
from the proton-nucleus elastic scattering data. In the calculation, we take the values of 1.76 MeV, 
500 c^/GeV^ for the Cmom and e, respectively, which result in the effective mass m* /m—0.75 in 
nuclear medium at saturation density for symmetric nuclear matter. The parameter x is changed 
as the strength of the mass splitting, and the values of -0.65 and 0.65 are respective to the cases of 
m* > m* and m* < m*. Shown in Fig. 1 is a comparison of different stiffness of nuclear symmetry 
energy after inclusion of the momentum-dependent interactions. One can see that all cases cross 
at saturation density with the value of 31.5 MeV. The local part of the symmetry energy can 
be adjusted to reflect the uncertain behavior of the symmetry energy at sub- and supra-normal 
densities. 
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The evolution of mesons (here mainly pions, kaons and etas) is also determined by the Hamil- 
tonian as follows: 

Hm =E(^.^°^^ + c^(p.,pO)- (9) 
i=i ^ ' 

Here the Coulomb interaction is given by 

^Coul ^ ^ ^IQ) 

i=i 

where the Nm and Nb are the total numbers of mesons and baryons including charged reso- 
nances. Here, we use the self-energy in vacuum for pions, namely c<;(pj,pj) = y^pf + with 
the momentum and the mass mu of the pions. The in- medium potentials for pion and kaon 
propagations can be chosen in the model [HI [T7]. We have included the resonances A(1232), 
A^*(1440), A^*(1535) in nucleon-nucleon (NN) collisions and their decaying into pions and etas. 
Furthermore, strange particles are directly created by inelastic hadron-hadron collisions [15]. The 
probability reaching a channel in a NN collision process is determined by its contribution of the 
channel cross section to the total cross section. The choice of the channel is done randomly by 
the weight of the probability. 

The preequilibrium nucleons in heavy-ion collisions are mostly produced during the compres- 
sion stage. Therefore, the high-density information of nuclear phase diagram can be extracted from 
the preequilibrium nucleon emissions, which are constrained with the transverse momentum and 
longitudinal rapidity distributions of free nucleons in this work. Shown in Fig. 2 is a comparison 
of the neutron/proton (n/p) yields of free nucleons within the mid-rapidity bin of \y/yproj\ <0.3 for 
the cases of hard and supersoft symmetry energies taken from the ^^^Sn+^^^Sn reaction at a beam 
energy of 400 MeV/nucleon with the mass splittings of m* > m* in the left panel and m* < m* 
in the right window, respectively. One can see that the distribution structure is sensitive to the 
collision centrality because of the different high-density reaction zone being formed, but weakly 
depends on the mass splitting. A flat spectrum at high transverse momentum is observed for the 
mass splitting of m* < m* because the more energetic neutrons are emitted in the case, which 
is consistent with the stochastic mean-field (SMF) calculations [18]. One expects to see a larger 
n/p value with a harder symmetry energy owing to a stronger repulsive (attractive) force acting 
on neutrons (protons) in the mean-field evolutions. It indeed exists in the near-central heavy-ion 
collisions. In particular, the difference of the hard and supersoft symmetry energies is larger at 
low transverse momentum. In the semi-central collisions, the influence of the symmetry energy 
exhibits an opposite distribution at high transverse momentum because of the secondary collisions 



with participant nucleons. And a flat spectrum appears for tlie supersoft symmetry energy. Tfie 
interest findings are different with the fermi-energy heavy colhsions [T2] . 

In order to reduce the systematic errors and the influence of Coulomb force, we calculated the 
transverse momentum distributions of the double ratios of free nucleons of reaction partners in 
^^^Sn+^^^Sn over ^^^Sn+^^^Sn collisions as shown in Fig. 3. One notices that the symmetry energy 
results in a variation of the value about 5% ~10%, in particular, at the low transverse momentum 
in near central collisions. The hard symmetry energy enhances the ratios because of the larger 
values at supra-saturation densities. The results are consistent with the IBUU04 calculations 
except for at high transverse momentum [12]. The mass splitting slightly affect the distributions. 
Experimental measurements are very anticipant to constrain the high-density symmetry energy 
from the preequilibrium nucleon emissions. Again, the influence of the symmetry energy on 
the double ratios in the case of peripheral collisions is opposite at high transverse momentum. 
Overall, the ratio increases with the transverse momentum in both collision centralities. Figure 
4 is a comparison of the double ^H/^He ratio of light clusters with hard and supersoft symmetry 
energies for the different mass splittings of m* > m* in the left window and m* < m* in the 
right window, respectively. The effect of the symmetry energy only appears at high transverse 
momentum in the semi-central collisions, but with larger error bars. 

Production of pions in heavy-ion collisions at near threshold energies is a primary product in 
nucleon-nucleon inelastic scattering and has been verified as a sensitive probe of nuclear symmetry 
energy [31 HI |5], which can be easily detected in experimentally and be mainly created in the 
domain at supra-saturation densities of compressed nuclear matter. The double n^/n^ ratios of 
total charged pions produced in the isotopic systems ^^^Sn-|-^^'^Sn over ^^^Sn-|-^^^Sn as a function of 
transverse momentum are calculated as shown in Fig. 5 at the beam energy of 400 MeV/nucleon 
for near-central collisions with the mass splittings of m* > m* (left panel) and m* < m* (right 
panel), respectively. One can see that the effect of symmetry energy is obvious at high transverse 
momentum. The hard symmetry energy leads to the larger values of the double ratios, in particular 
at high transverse momentum, which results from the fact that a number of tt^ are produced with 
the hard symmetry energy and pronounced in neutron-rich systems in the model O [H]. The 
conclusions does not satisfy the results in Ref. [12], where a softer symmetry energy leads to a 
larger value of the double ratio for the energetic pions and an opposite distribution appears at 
low kinetic energies. The pion production is complicated not only involving the nucleon-nucleon 
inelastic collisions such as the A^A^ — )■ A^A and the inverse channel A^A — )■ A^A^, but also being 
related to the absorbtion by a nucleon. The mass splitting slightly influences the high-momentum 



distribution, but does not change the effect of symmetry energy. 

In summary, within the transport model LQMD we have investigated the transverse emis- 
sions of free nucleons, Ught chistcrs and charged pions produced in isotopic reaction systems of 
^^^Sn+^^^Sn and ^^^Sn+^^^Sn at the beam energy of 400 MeV/nucleon. The single and dou- 
ble ratios of isospin particles are analyzed with different collision centralities and effective mass 
splittings. It is found that the single and double neutron/proton ratios of mid-rapidity nucleon 
emissions and the double 7r~/7r"'" yields are sensitive to the density dependence of the symmetry 
energy, which are useful probes for constraining the high-density behavior of the nuclear symmetry 
energy. The collision centrality changes the high- momentum distributions of the single and double 
ratios, in particular, for the free nucleon emissions. However, the mass splitting slightly affects 
the transverse emission of the isospin ratios. 
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Figure 1: Nuclear symmetry energy as a function of baryon density with the MDI interaction for 
the different stiffness (supersoft, soft, hnear and hard). 
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Figure 2: Transverse momentum distributions of neutron/proton ratio of free nucleon emis- 
sions within the rapidity bin \y/yproj\ <0.3 in the ^^^Sn+^^^Sn reaction at incident energy of 
400 MeV/nucleon for near-central (b=l fm) and semi-central (b=5 fm) collisions with the mass 
splittings of m* > m* (left window) and m* < m* (right window), respectively. 



in 




Figure 3: The double neutron/proton ratio of free nucleons taken from the reactions of ^^^Sn+^^^Sn 
and ^^'^Sn+^^^Sn at an incident energy of 400 MeV/nucleon for different mass sphttings of m* > m* 
(left window) and m* < m* (right window), respectively. 
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Figure 4: The same as in Fig. 3, but for the triton and hehum-3 emissions. 
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Figure 5: Comparison of the transverse momentum dependence of the double 7r~/7r+ ratio for 
near-central (b=l fm) ^^^Sn+^^^Sn and ^^^Sn+^'-^^Sn collisions at the incident energy of 400 
MeV/nucleon with the hard and supersoft symmetry energy, but for different mass splittings. 



